BEST AVAI LASlfLl! ( 0©$Y 



to 

CM 
CO 



Q. 

LU 



(19) 



J) 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(12) 



(11) EP 1 324 455 A1 

EUROPEAN PATENT APPLICATION 



(43) Date of publication: 

02.07.2003 Bulletin 2003/27 

(21) Application number: 01811253.2 

(22) Date of filing: 21 .12.2001 



(51) mtci7: H 02 H 7/26, H02H 3/40, 
H02J 3/24 



(84) Designated Contracting States: 


(72) Inventors: 


AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 


♦ ZIma, Marek 


MC NL PT SE TR 


5400 Baden (CH) 


Designated Extension States: 


• Rehtanz, Christian 


AL LT LV MK RO SI 


5405 Baden (CH) 


(71) Applicant: ABB RESEARCH LTD. 


(74) Representative: ABB Patent Attorneys 


8050 Zurich (CH) 


c/o ABB Schweiz AG 




Brown Boveri Strasse 6 




5400 Baden (CH) 



(54) Electric power transmission network state estimation 



(57) In a method, computer program and system for 
estimating a state of an electric power transmission net- 
work, time-stamped phasor measurements are made by 
phasor measurement units (1) located at a first set of 
locations (11,12) throughout the network, and a plurality 
of voltage and/or current phasors at a second set of lo- 
cations is estimated, where the locations of the first 
(11,12) and second sets of locations are distant from 
one another, and a status of network elements 
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(L4,L5,13) that are remotely located from the phasor 
measurement units is inferred from the phasor meas- 
urements and the estimated phasors. 

In a preferred variant of the invention, the network 
elements are transmission lines (L1...L8), switchgear, 
under load tap changing transformers, compensation 
devices (1 3) or FACTS devices, and the status informa- 
tion describes whether a network element (L4,L5,13) is 
electrically connected to or disconnected from the net- 
work. 
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Description 

Field of the Invention 

[0001] The invention relates to large-scale electric 5 
power transmission networks, and, more particularly, to 
a method, a computer program and a system for esti- 
mating a state of an electric power transmission network 
according to the preamble of claims 1 , 8 and 9, respec- 
tively. 10 

Background of the Invention 

[0002] Electric power transmission and distribution 
systems or networks comprise high-voltage tie lines for 15 
connecting geographically separated regions, and sub- 
stations for transforming voltages and for switching con- 
nections between lines. Power generation and load flow 
in a network with several substations is controlled by an 
energy management system. For managing the net- 20 
work, it is desirable to determine a state of the network, 
in particular load flows and stability margins. Measure- 
ments of average RMS (root mean square) values of 
voltages, currents, active power and reactive power 
flowing in the network and/or measurements of voltage 25 
and current phasors are made at various places in the 
network and collected by substation automation (SA) 
systems and by a system control and data acquisition 
(SCADA) system. The SA and SCADA systems further 
provide information about a status of network elements 30 
such as switchgear, compensation devices and variable 
transformers. The status of switchgear (i.e. circuit 
breakers, disconnectors, etc.) allows to infer the topol- 
ogy of the network, i.e. how the lines and devices are 
connected, and which lines, power generators and con- 35 
sumers are disconnected. Additionally, measurements 
of voltage and current phasors are made and collected. 
Ail the abovementioned measurements from throughout 
the network are collected at the energy management 
system, providing a snapshot of the state of thenetwork. *o 
The snapshot allows to compute the state of the net- 
work, as described e.g. in "State estimation with phasor 
measurements", Phadke et. al., IEEE Transactions on 
power systems, Vol. PWRS-1 , No. 1 , Feb. 1 996. There, 
all positive sequence bus voltages and some positive 45 
sequence currents are measured, the network topology 
is assumed to be known, and unknown voltages, phase 
angles and power flows are computed. 

Description of the Invention 50 

[0003] It is an object of the invention to create a meth- 
od, a computer program and a system for estimating a 
state of an electric power transmission network of the 
type mentioned initially, which uses less measurements, 
[0004] These objects are achieved by a method, a 
computer program and a system for estimating a state 
of an electric power transmission network according to 



the claims 1 , 8 and 9, respectively. 
[0005] In the inventive method, computer program 
and system for estimating a state of an electric power 
transmission network, time-stamped phasor measure- 
ments are made by phasor measurement units located 
at a first set of locations throughout the network, and a 
plurality of voltage and/or current phasors at a second * 
set of locations is estimated, wherein the locations of ' 
the first and second sets of locations are distant from 
one another, and a status of network elements that are 
remotely located from the phasor measurement units is 
inferred from the phasor measurements and from the 
estimated phasors. 

[0006] This allows to estimate a network state without 
having to collect status information of network elements 
such as switchgear, compensation and FACTS devices 
and variable transformers through a SCADA or a Sub- 
station Automation system. 

[0007] An additional advantage is conferred by the 
fact that phasor measurements incorporate high reso- • J 
lution timing information. It therefore is also possible to 
obtain a much better temporal resolution for the state 
estimation, since conventional SCADA systems do not 
provide time stamped information about e.g. switchgear 
status. This better temporal resolution allows faster and 
better quality control of the network. 
[0008] In a preferred variant of the invention, the net- 
work elements are transmission lines, switchgear, under 
load tap changing transformers, compensation devices 
or FACTS devices, and the status information describes 
whether a network element is electrically connected to 
or disconnected from the network and optionally also 
describes element specific parameters such as a trans- 
former tapping position or a FACTS impedance value. 
For switchgear such as breakers or disconnectors, the 
status information describes whether a switch is open 
or closed. 

[0009] In a further preferred variant of the invention, 
the state that is estimated comprises current and load ' > 
flows, and optionally indicators characterising network 
stability. This is possible, since the invention allows to 
determine, from phasor measurements alone, which 
network elements are connected, i.e. the network topol- 
ogy. No separate measurements or status information 
about the network topology are required. Therefore, the 
invention can work completely independent of and re- 
dundant to existing SCADA systems. 
[001 0] The system for estimating a state of an electric 
power transmission network according to the invention 
comprises a data processing device configured to com- w 
pute a plurality of estimated voltage and/or current 
phasors at a second set of locations, the locations of the 
first (11 ,12) and second sets of locations being distant < 
from one another, and to compute, from the phasor 
measurements and from the estimated phasors, a sta- 
tus of network elements (L4,L5,1 3) that are remotely lo- 
cated from the phasor measurement units. 
[0011] The computer program for estimating a state 
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of an electric power transmission network according to 
the invention is loadable into an internal memory of a 
digital computer, and comprises computer program 
code means to make, when said computer program 
code means is loaded in the computer, the computer ex- 
ecute the method according to the invention. In a pre- 
ferred embodiment of the invention, a computer pro- 
gram product comprises a computer readable medium, 
having the computer program code means recorded 
thereon. 

[0012] Further preferred embodiments are evident 
from the dependent patent claims. 

Brief Description of the Drawings 

[0013] The subject matter of the invention will be ex- 
plained in more detail in the following text with reference 
to preferred exemplary embodiments which are illustrat- 
ed in the attached drawings, in which: 

Figure 1 schematically shows part of an electric pow- 
er system, 

Figure 2 shows an equivalent circuit for a transmis- 
sion line 

Figure 3 schematically shows an electrical substa- 
tion 

Figure 4 shows an equivalent circuit for a tap chang- 
ing transformer 

Figure 5 schematically shows the structure of a sys- 
tem according to the invention. 

[0014] The reference symbols used in the drawings, 
and their meanings, are listed in summary form in the 
list of reference symbols. In principle, identical parts are 
provided with the same reference symbols in the figures . 

Detailed Description of Preferred Embodiments 

[0015] Figure 1 schematically shows part of an high- 
or medium voltage electric power system or network 
comprising nodes n1 ... n5 such as busbars, and trans- 
mission lines L1 ... L.6. A plurality of nodes n1 , n2 is lo- 
cated at a first substation 11 , another node n5 is located 
at a second substation 12. A compensation element 13 
is connected to a fourth node n4. At the first and second 
substation 11,12, phasor measurement units (PMU) 1 
are configured to measure voltage and/or current 
phasors at the lines L1 ,L2,L3,L6 leaving the respective 
substation 11,12. A third node n3 and the fourth node 
n4 are not equipped with PMU's. The third node n3 and 
fourth node n4 are located remote or geographically dis- 
tant from one another and from the first substation 1 1 
and the second substation 12. They are separated by 
the lengths of the lines L1 ...L6 connecting them, typical- 
ly at distances between one kilometre and several tens 
or hundreds of kilometres. Where the lines meet the 
nodes, switchgear devices, not shown in the figures, are 
configured for electrically disconnecting lines. 



[0016] Phasor data is determined with a phasor 
measurement unit residing, for example, at a feeder at 
the bay level of substations or at branching points along 
transmission lines. A voltage phasor represents, for ex- 
5 ample, a voltage of the feeder or line, while a current 
phasor represents current flowing through the feeder or 
line. 

[001 7] The phasor data represents a phasor and may 
be a polar number, the absolute value of which corre- 

10 sponds to either the real magnitude or the RMS value 
of a quantity, and the phase argument to the phase an- 
gle at zero time. Alternatively, the phasor may be a com- 
plex number having real and imaginary parts or the 
phasor may use rectangular or exponential notation. By 

is contrast, conventional sensing devices used in power 
networks generally measure only scalar, average rep- 
resentations, such as the RMS value of a voltage, cur- 
rent etc. 

[0018] The phasor data is collected from phasor 

20 measurement units that are distributed over a large ge- 
ographic area, i.e. over tens to hundreds of kilometres. 
Since the phasor data from these disparate sources are 
analysed in conjunction, they must refer to a common 
phase reference. Therefore, the different phasor meas- 

25 urement units must have local clocks that are synchro- 
nised to within a given precision. Such a synchronisation 
of the phasor measurement units is preferably achieved 
with a known time distribution system, for example the 
global positioning (GPS) system. In atypical implemen- 

30 tation, the phasor data 9 is determined at least every 
200 or every 100 or preferably every 40 milliseconds, 
with a temporal resolution of preferably less than 1 mil- 
lisecond. In a preferred embodiment of the invention, the 
temporal resolution is less than 1 0 microseconds, which 

35 corresponds to a phase error of 0.2 degrees. Each 
measurement is associated with a time stamp derived 
from the synchronised local clock. The phasor data 
therefore comprises time stamp data. 
[0019] According to the invention, it is possible to de- 

^o termine a status of the fourth line L4 and the fifth line 
L5, i.e. whether said lines are connected or disconnect- 
ed, without measurements or SCADA data about a sta- 
tus of switchgear at the third node n3 or the fourth node 
fourth node n4, and even without phasor measurements 

45 made at the third or fourth node n3,n4. The term "status" 
is used to denote properties of network elements such 
as "connected/disconnected", "open/closed", "tap ratio 
value" that usually are set by a control device and an 
actuator associated with a network element, 

so [0020] The term "state" comprises both the "status" of 
network elements as well as e.g. current and power 
flows throughout the network, and optionally closeness 
to instability. The state is a result of electrical parameters 
of the network and of network topology, which is deter- 

55 mined by the status of network elements. 

[0021] Figure 2 shows a standard it-equivalent circuit 
for a transmission line or line segment. V1 is a voltage 
phasor at the first node n1 , 111 is a phasor of a current 
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flowing out of the first node n1 into the first line L1 . Both 
V1 and 111 are measured by a PMU 1 located at the first 
node n1 . The first line L1 is represented by a line admit- 
tance Y1 and shunt admittances Y11 and Y31, with 
Y11=Y31. By convention, for a current, the first index 
denotes the node from which the current flows, and the 
second index denotes the line into which it flows. For a 
shunt admittance, the first index denotes the node at the 
side of the line where the shunt lies, and the second in- 
dex denotes the line. 

[0022] V3 is an unknown voltage phasor at the third 
node n3. V3 is calculated as 

» n (Y1+Y11)V1-111 

[0023] In a similar fashion, referring to the power sys- 
tem shown in Fig. 1 , the voltage V4 at the fourth node 
n4 is calculated as 
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15 



V4 = 



(Y6+Y56)V5-l56 
Y6 



from the voltage and current phasors V5,I56 measured 
at the f ifth node n5. 

[0024] The currents flowing through the first, second 
and third line L1 ,L2,L3 into the third node n3 are calcu- 
lated from V3 and the measurements. For example, 

I31=(Y1+Y31)V3-Y1 V1 

[0025] The total current from the third node n3 to the 
fourth node n4 therefore is 



termine which of the two is disconnected. This is, how- 
ever, not relevant for a load flow calculation that uses 
this topology information - exactly because the line ad- 
mittances are equal, regardless of which of the two lines 
L4,L5 is disconnected. 

[0027] In order to determine a status of the compen- 
sation element 13 that is not observed directly, i.e. it is 
not equipped with a PMU 1 , the procedure for determin- 
ing voltages and detecting line status as shown above 
is executed. Then, all currents flowing into the fourth 
node n4 are computed, depending on the network to- 
pology. For example, if both the fourth and fifth line L4,L5 
are connected, and from 

I45 = (Y5 + Y45) V4 - Y5 V5 

it follows that the sum of currents flowing into the fourth 
node n4 through lines L4,L5 and L6 is 



20 



25 



30 



35 



l3to4 = - ( l31 + 132 + ,33 )' 

If this current is zero, then the status of both the fourth 
and fifth line L4,L5 is "disconnected". If it is nonzero, it 
is necessary to determine which of said lines L4.L5, if 
any, is disconnected. Assuming that only the fourth line 
L4 is connected and that the fifth line L5 is out of service, 
the current flowing into the fourth line L4 from the third 
node n3 can be calculated as 

134 = (Y4 + Y34) V3 - Y4 V4 

which ought to be equal to l 3to4 . If it is equal, then the 
assumption is valid. It is not equal, then the inverse as- 
sumption "only the fifth line L5 is connected" is made 
and checked. If this assumption is invalid as well, then 
both lines L4,L5 are connected. 

[0026] In the unlikely case that the admittances of the 
fourth line L4 and the fifth line L5 are be equal, it is only 
possible to determine whether none, one or both of 
these lines L4.L5 is connected. It is not possible to de- 
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45 



I 4in = -(I45 + I44+ 146). 

If l 4in is zero, then the compensation element 13 is dis- 
connected. If l 4in is not zero, the compensation element 
13 is in use. 

[0028] For FACTS (flexible AC transmission) devices, 
with knowledge about the internal workings of the de- 
vice, status information may be determined from the cur- 
rent and voltage phasors at the device. Regardless of 
the internal workings, whether a FACTS device is con- 
nected or not is inferred in the same way as shown 
above for the compensation element 1 3. 
[0029] Figure 3 schematically shows an third electri- 
cal substation 13 with busbars acting as nodes n6,n7, 
a transformer 31 connected to a sixth node n6, a busbar 
connector 32 connecting the sixth node n6 to a seventh 
node n7, and lines L7.L8 and a load 33, each line and 
load 33 being connected to one of the nodes n6,n7 and 
having an associated PMU 1 . Although the busbar con- 
nector 32 is located in the third substation 13, no direct 
observation with a PMU 1 or current transducer is avail- 
able. Therefore, the status of the busbar connector 32 
is inferred from a currents balance, using phasor meas- 
urements at eitherthe sixth node n6 orthe seventh node 
n7. For example, if 
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55 



tirrtoe = " ( ^ansformer + ,67 ) 

is zero, then the busbar connector 32 is open. Other- 
wise, it is closed. 

[0030] Transformers without automatic tap changing 
are modelled like static line elements, i.e. by symmetri- 
cal rc-equivalent circuits, and have only an on/off status. 
[0031] Figure 4 shows an equivalent circuit for an un- 
der load tap changing (ULTC) transformer arranged be- 
tween nodes with voltage phasors V1 and V2 and a cur- 
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rent phasor 12 flowing into the transformer from the node 
with vottage V2. The admittances of the circuit depend 
on a tap ratio t. The ULTC transformer is equipped with 
a control device that adjusts the voitage ratio of the 
transformer within certain limits around a nominal volt- 
age ratio. The control device tries to drive a regulated 
voltage, in this case V2, at one side of the transformer 
to within a deadband around a given nominal value 
V2nom. This adjustment is represented by the tap ratio 
t. For t=1 , the transformer is set to the nominal voltage 
ratio. Typically, the tap and voltage ratio can be varied 
by 5% or 10% around the nominal value. In order to de- 
termine the current and load flows of the power network, 
the tap ratio t must be known. According to the invention, 
it is not necessary to obtain information from the local 
ULTC transformer's control device. Instead, the tap ratio 
t is deduced from phasor information. The phasor infor- 
mation is determined either by one or more PMUs 1 lo- 
cated close to the ULTC transformer, or is inferred from 
remotely located PMUs 1 as described above. Likewise, 
it is determined whether the ULTC transformer is con- 
nected at all. Parameters of the ULTC transformer that 
must be known are its nominal admittance Y, incremen- 
tal tap ratio, minimal and maximal tap ratio t and a dead- 
band area around the nominal regulated value. 
[0032] Two cases are possible: Either the actual or 
present value V2nom of the regulated voltage is known, 
or it is not known. 

[0033] For the first case, the tap ratio t is determined 
from the regulated voltage V2 alone. Let I. ..I denote the 
absolute value of a complex variable. Then, 

if IV2l<0.9IV2noml => t = 0.9 



if IV2l>1.1IV2nom1 =>t= 1.1 



if 0 98 < ,wo . < 1 .05 => t = 1 
IV2noml 

where 0.98 and PMU 1 .05 are bounds of the deadband, 
and in all other cases, the tap ratio t is rounded to the 
nearest possible value, according to the incremental tap 
ratio. 

[0034] The admittances of the equivalent circuit are 
calculated from their nominal values and from the tap 
ratio t as 

y ultc = x Y 
Y 1ULTC =t(M)Y 

Y 2ULTC = 0-t)Y 
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[0035] For the second case, when the nominal value 
V2nom of the regulated voltage is not known, the volt- 
age phasors V1 ,V2 on both sides of the ULTC trans- 
former, and the current phasor 12 on the side of the reg- 
s ulated voltage V2 must be known. The tap ratio t is cal- 
culated as 

_ Y V2-I2 
YV1 

10 

[0036] For a given network topology, the placement 
of PMUs 1 and the exact sequence of operations for 
computing voltages and inferring the status of network 
element status is determined e.g. manually and with 
* 5 computer assistance in an engineering phase of the sys- 
tem according to the invention. Preferably, a set of 
PMUs 1 that is required for other purposes is given. 
Starting with their measurements, equations for deter- 
mining as many voltages and current phasors at nodes 
and lines not equipped with PMUs are generated. This 
is done for all possible topologies, i.e. for all network 
elements being either connected or disconnected. If, for 
certain cases, the existing phasor measurements are 
not sufficient, an additional PMU 1 is introduced either 
at the location where information is required, or at a 
node connected to said location. 
[0037] In a preferred embodiment of the invention, the 
placement of PMUs 1 is minimised in a two-step proce- 
dure: in a first step, the minimum number of substations 
equipped with PMUs 1 is determined, assuming that in 
such a substation, each node equipped with an instru- 
ment transformer has an associated PMU. In a second 
step, a minimum number of PMU's for each substation 
is determined. 

[0038] Figure 5 schematically shows the structure of 
a system according to the invention. Each of a plurality 
of PMUs 1 transmits current and/or voltage phasor in- 
formation to a network element status detection 2. The 
network element status detection 2 is a data processing 
device configured to determine the status of network el- 
ements, as shown above. It is configured to receive, as 
an additional input, network element parameters 3 such 
as line and shunt admittances, tap ratio increments, lim- 
its and deadband values, regulated voltages etc. In a 
preferred embodiment of the invention the network ele- 
ment status detection 2 transmits the status of network 
elements together with the phasor information as status 
and phasor information 4 to a state calculation/estima- 
tion data processing device 5. This information 4 is used 
by the state calculation/estimation 5 to estimate a state 
of the network, e.g. by doing a load flow analysis, in 
which topology information is given by the status of net- 
work elements without the need of additional measure- 
ments. 
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List of designations 


[0039] 




1 


phasor measurement unit (PMU) 


2 


network element status detection 


3 


network element parameters 


4 


status and phasor information 


5 


state calculation/estimation 


11 


first substation 


12 


second substation 


13 


compensation element 


13 


third substation 


31 


transformer 


32 


busbar connector 


33 


load 


n1, n2...n7 


first node, second node, ... seventh 




node 


L1, L2, ... L8 


first line, second line, ... eighth line 




current phasor 


v.. 


voltage phasor 


Y.. 


complex admittance 


t 


tap ratio 



25 

Claims 

1 . A method for estimating a state of an electric power 
transmission network, in which time-stamped 
phasor measurements are made by phasor meas- 30 
urement units (1) located at a first set of locations 
(11,12) throughout the network, characterized in 
that 

a plurality of voltage and/or current phasors at a 
second set of locations is estimated, and wherein 35 
the locations of the first (11 ,12) and second sets of 
locations are distant from one another, and 
that a status of network elements (L4,L5,13) that 
are remotely located from the phasor measurement 
units (1 ) is inferred from the phasor measurements 40 
and the estimated phasors. 

2. The method according to claim 1 , in which voltages 
and/or currents at a substation that is not equipped 
with phasor measurement units (1) is estimated. 45 

3. The method according to one of the preceding 
claims, in which a state of the network comprising 
load and current flows of the network is estimated 

f rom the measured and estimated voltages and cur- so 
rents and from an estimated topology of the net- 
work. 

4. The method according to one of the preceding 
claims, in which the network elements comprise 55 
switchgear and the status of network elements 
comprises information on whether a switch is open 

or closed. 



5. The method according to one of the preceding 
claims, in which the network elements comprise 
transmission lines (L1...L8) and the status of net- 
work elements comprises information on whether a 
transmission line is disconnected or in service. 

6. The method according to one of the preceding 
claims, in which the network elements comprise un- 
der load tap changing transformers and the status 
of network elements comprises information on a tap 
ratio of the transformer. 

7. The method according to one of the preceding 
claims, in which the network elements comprise at 
least one compensation device (13) or FACTS de- 
vice and the status of network elements comprises 
information on whether said device is electrically 
connected to or disconnected from the network. 

8. A computer program for estimating a state of an 
electric power transmission network loadable into 
an internal memory of a digital computer, compris- 
ing computer program code means to make, when 
said program is loaded in the computer, the compu- 
ter execute a procedure for processing time- 
stamped phasor measurements originating from 
phasor measurement units (1) located at a first set 
of locations (11 ,12) throughout the network, 
characterized in that the program is configured to 
compute a plurality of estimated voltage and/or cur- 
rent phasors at a second set of locations, the loca- 
tions of the first (11,12) and second sets of locations 
being distant from one another, and to 
compute, from the phasor measurements and from 
the estimated phasors, a status of network ele- 
ments (L4,L5, 1 3) that are remotely located from the 
phasor measurement units. 

9. A system for estimating a state of an electric power 
transmission network from time-stamped phasor 
measurements that originate from phasor measure- 
ment units (1) located at a first set of locations 
(11,12) throughout the network, characterized In 
that 

the system comprises a data processing device 
configured to compute a plurality of estimated volt- 
age and/or current phasors at a second set of loca- 
tions, the locations of the first (11,12) and second 
sets of locations being distant from one another, 
and to 

compute, from the phasor measurements and from 
the estimated phasors, a status of network ele- 
ments (L4.L5, 1 3) that are remotely located from the 
phasor measurement units. 

10. The system according to claim 9, wherein the data 
processing device is configured to estimate voltag- 
es and/or currents at a substation that is not 
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equipped with phasor measurement units (1). 
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